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(54) Signal bias removal for robust telephone speech recognition. 

(57) A signal bias removal (SBR) method based on the maximum likelihood estimation of the bias for 
minimizing undesirable effects in speech recognition systems is described. The technique is readily 
applicable in various architectures including discrete (vector-quantization based), semicontinuous and 
continuous-density Hidden Markov Model (HMM) systems. For example, the SBR method can be 
integrated into a discrete density HMM and applied to telephone speech recognition where the 
contamination due to extraneous signal components is unknown. To enable real-time implementation, 
a sequential method for the estimation of the bias (SSBR) is disclosed. 
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Technical Field 

This invention relates to a method for removing a signal bias from a transmitted speech signal. In particular, 
the present invention is a signal bias removal method which may be integrated into a Hidden Markov Model 
5 (HMM)-based speech recognition system to minimize the effects of unknown adverse conditions that typically 
contaminate the speech signal in a telephone channel. 

Background of the Invention 

10 A speech signal transmitted through a telephone channel often encounters unknown variable conditions 

which significantly deteriorate the performance of state-of-the-art Hidden Markov Model (HMM)- based speech 
recognition systems. Undesirable components due to ambient noise and channel interference, as well as dif- 
ferent sound pick-up equipment and articulatory effects, render such recognition systems unsuitable for many 
real- world applications. 

15 Noise is usually considered to be additive to the speech signal. The spectrum of a real noise signal, such 

as that produced from fans and motors, is generally not flat and can often cause a considerable degradation 
in the performance of a speech recognizer. 

Channel interference, both linear and non-linear, can also have a serious impact on a speech recognizer. 
An effect of a typical telephone channel is that it band pass filters the transmitted signal between 200 Hz and 

20 3200 Hz, with variable attenuations across the different spectral bands. If this filtering action is not made con- 
sistent when training and testing a speech recognizer, severe consequences on the performance may result, 
tn addition, the use of different microphone transducers can create an acoustic mismatch in the training and 
the testing conditions. 

Another source of degradation in the performance of a speech recognizer pertains to articulation effects. 

25 Changes in articulation usually occur due to environmental influences (known as the "Lombard effect"), but 
may occur merely when speaking to a machine. Articulation effects are a major concern in telephone speech 
recognition, especially in situations where, for example, a customer is talking to an automatic speech recog- 
nizer from a public phone-booth situated near a major highway. 

Prior art efforts to minimize extraneous signal components for robust speech recognition have centered 

30 upon three major areas. First, processing the speech signal to remove an estimate of the noise. Typical ex- 
amples include spectral subtraction, cepstral normalization, noise masking and robust feature analysis. Sec- 
ond, adapting the recognizer's models to noise without modifying the speech signal. Third, applying a robust 
distortion measure that emphasizes the regions of the spectrum that are less corrupted by noise. 

35 Summary of the Invention 

The present invention addresses the problem of unknown signal bias contamination of a speech signal by 
introducing a signal bias removal (SBR) method into a speech recognition system. The SBR method separates 
an estimate of the bias signal from the input speech signal. This differs from prior art methods that are designed 

40 to normalize or to correct the cepstral vectors. In addition, the separation formulation can be done in the ob- 
servation space as welt as in the probability space. 

The SBR method may be utilized during both training and testing in a speech recognition system. After 
performing feature analysis on a training speech signal, an estimate of the bias is computed based on maxi- 
mizing a likelihood function. Next, the estimate of the bias is subtracted from the speech signal to arrive at a 

45 tentative speech value. Computing the estimate of the bias and subtracting the estimate from the signal are 
repeated a predetermined number of times, and each iteration uses the previous tentative speech value to com- 
pute the next bias estimate to arrive at a reduced bias speech signal value. Next, a codebook of centroids is 
generated, and then the estimate of the bias is recomputed and subtracted from the tentative speech signal 
again until an optimal set of centroids are generated. The reduced bias speech signal value and optimal cen- 

so troids are then used as the training input to the speech recognizer. The testing phase then consists of utilizing 
the codebook generated during training to compute an estimate of the bias for each utterance based on max- 
imizing a likelihood function, subtracting the estimate of the bias from the speech signal to obtain a tentative 
speech value, and then repeating these two steps a preset number of times to result in a reduced bias speech 
value. The reduced bias speech value is then used as the input for the speech recognizer. 

ss A further aspect of the present invention is that the SBR method may be applied after a speech recognition 

system has undergone any suitable training phase. Thus, SBR can be used in the testing phase only, utilizing 
the set of centroids generated during the training phase. 

Yet another aspect of the invention concerns a sequential SBR method (SSBR) which permits real-time 
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implementation of the invention on current platform recognizers without any major structural change. The 
SSBR method enables the processing of the bias at the frame level, rather than at the utterance level, without 
imposing a look-ahead frame delay. 

The S8R and SSBR methods are readily applicable to various HMM architectures, such as discrete (Vector 
5 Quantization-based), semi-continuous and continuous density HMM systems. In addition, for both the SBR 
and the SSBR methods, bias removal is carried out as an independent process following feature analysis and 
preceding recognition. Thus, the present invention may be integrated into a discrete density HMM system used 
for telephone speech recognition. 

10 Brief Description of the Drawings 

Fig. 1 is a block diagram of a distorted telephone network; 

Fig. 2 A is a block diagram illustrating the integration of the signal bias removal technique of the present 
invention in a speech recognition system; 
15 Fig. 2B is a block diagram illustrating the integration of the signal bias removal technique of the present 

invention in an HMM-based recognition system; 

Fig. 3 is a flowchart illustrating the use of the invention when training a speech recognition system; 
Fig. 4 is a flowchart illustrating the use of the invention when testing an input speech signal In a speech 
recognition system; 

20 Fig. 5A is a block diagram illustrating an apparatus incorporating the present invention; 

Fig. 5B is a block diagram illustrating another apparatus incorporating the present invention; 

Fig. 6 is a table illustrating the percentage word error and rates of insertion, deletion and subtraction for 

a baseline system with either SBR or CMS; 

Fig. 7 is a plot of the norm of the cepstral bias averaged over the training and the test data at every iteration 
25 of the SBR method; 

Fig. 8 is a plot of the quantization error averaged over the training and the test data at every iteration of 
the SBR method; 

Fig. 9 is a histogram for the second bias coefficient when testing on a first database (DB1) ; 
Fig. 1 0 is a histogram for the second bias coefficient when testing on a second database (DB2) ; 
30 Fig. 11 is a histogram of the second bias coefficient when testing on DB2 following ten iterations of the 

SBR method; 

Fig. 12 is a table illustrating the percentage word error for the baseline system with SBR, SSBR, or CMS; 
Fig. 1 3 is a plot of the word error rate as a function of the codebook size when using SBR; 
Fig. 14 is a plot of the word error rate as a function of the codebook size when using SSBR; and 
35 Fig. 1 5 is a table illustrating the percentage word error for different string lengths before and after utilizing 

SBR. 

Detailed Description 

40 Fig. 1 is a schematic block diagram of a distorted telephone network 1 . A telephone speech signal X(a>) 

encounters a distortion effect having a multiplicative component, H(a>), due to distortion in the telephone chan- 
nel 2, and an additive component, N(a>), representative of the ambient noise. If X(a>) is the power spectrum of 
the original speech signal, then the received contaminated signal, Y(co), is modeled as: 

V(©) - H(co) + /V(g>), 

45 where H(<o) and N(a>) are "biases" which are assumed to be relatively constant throughout each utterance se- 
quence. The present invention estimates these biases and minimizes, or removes, their effects from the con- 
taminated signal Y(co). 

Fig. 2Ais a schematic block diagram of a speech recognition system 20 incorporating the present invention. 
The distorted signal Y((d) from a telephone channel is input to a feature analysis block 22, which performs a 

so sequence of feature measurements to form a "test pattern". The feature measurements are typically the output 
of any of several known spectral analysis techniques, such as filter bank analysis or a linear predictive coding 
(LPC) analysis. Typically, a bank of bandpass filters is used to screen the speech signal, and then a micropro- 
cessor is used to process the filtered signals. The results of the feature analysis are a series of vectors that 
are characteristic of the time-varying spectral characteristics of the speech signal. A codebook of these distinct 

55 analysis vectors is usually generated by one or more microprocessors utilizing a vector quantization (VQ) tech- 
nique. Vector quantization is known in the art, and is sometimes used as a preprocessor step to perform pre- 
liminary recognition decisions in order to reduce the computational load of a recognizer 24. 

Referring again to Fig. 2A, the output signal of the feature analysis block 22, V(<o), is used as the input to 
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P7 Th J ^ TToif BR> WOCk 26> The SBR blOCk 26 COmprises a codebook 25 and a bias removal section 
t.LT.H ®. BR f bl ° ck 26 f lrst «»««>•*»■ an Climate of the bias in the speech signal, and then removes the es- 
L?cn Z f 9e " erate 4 a si9nal for in P ut to • speech recognizer 24. The recognizer 24 characterizes 
he spectral properties of the frames of the speech pattern to output a text signal X"(o>)), which approximates 
the ongma I speech signal X(») with the bias removed. One well-known and widely used recognition technique 
mSX™ ^ °" e W 010,8 micro P roc «»*H>rs to perform speech recognition is the Hidden Markov Model 

S^f/T"? °' 8 eneratin 9 the text out P"t- present invention provides the recognizer 24 with a speech 
signal Y (<d) which contains a reduced amount of noise, thus permitting its output X'(<b)) to more closely ap- 
proximate the input speech signal X(o>). «oseiy ap 
Fig. 2B illustrates the integration of the SBR method 28 in a discrete density HMM recognizer system 21 

w" °T ' EH?" 8 renl0Val iS ° arried ° Ut aS a " inde P ende " t P^ess following feature analysis 23 and pre- 
nSS^* recognizer 29. This type of architecture enables the SBR technique to become an integrated 
part of the discrete density HMM. and the same method can be extended to or used as a front end processor 
for other HMM structures. 

The following derivation of the SBR method considers only the multiplicative spectral bias H(<o) however 
he process of dealing with both biases. H(co) and N(<o), relies on the same algorithmic development. Further' 
he formulation of the SBR method can be applied utilizing both spectral analysis and cepstral analysis. Ceps^ 
tral coefficients are the coefficients of the Fourier transform representation of the log magnitude spectrum 
and are known to provide a reliable feature set for speech recognition. 

The SBR method is based on maximizing the likelihood of a speech model in which the bias is considered 
as the unknown parameter. The likelihood function is defined as: 



where 

is an observation sequence of T frames. The speech modVi's^ 1-1.2. ... M ), where X, is the Markov 
model for a speech un.t V. The index t denotes the frame number ranging from 1 to T. Note that if the Markov 
chain in A is of zero* order or equiprobable, then the likelihood function defined above is equivalent to a nearest 

" 7 Q1 m ° del - Furthermore ' in its sim P |est ^ * can be just a codeword 
or centroid in the VQ codebook, where o> is assumed to be "bias free". 

Assuming an additive bias term, b, then with a simplified notation: 

y t = x t + b 

and 

... v={yi.y 2 ....y„ ...y T h 

and then 

P(V|d) = p(Y - b). 

The likelihood function thus becomes: 



p(y|^A)=n^ p(yt -^ i); 

and the maximum likelihood bias estimator, b. is the one that achieves: 



p[y\b.h) ^p(y\ b , A ) . 

The solution for the maximum likelihood bias estimate, b, can be found by using an iterative procedure 
consider a Gaussian local observation: 

Pfy t \b,X,) = K)exp{- |{(y ( - b) - n ( K£ 7 1 [(y, - b) - nj}, 
where X, = („,. z.) = (mean, covariance), and is the normalizing constant which does not depend on the bias 
b. Note that independence among the features is assumed, and thus Z, =l, the identity matrix. With an existing 
bias vector b, each adjusted observation is: 
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c = y* - t> 



and its nearest neighbor can be solved for, such that: 

5 max 

z t -\ii = argr m p[y e \b, A,) 
J J 

max 

J 

The likelihood function thus becomes: 

15 

p(Y\b,A) =*rexp{--t£ (y t -b-z t )>-{y t - b -z t )\ i 

20 where K is a constant that does not depend on the bias, b. By maximizing the quadratic function for p(Y \ bj\) 
with respect to the bias vector b, a unique solution for the bias estimate is guaranteed: 



J t-1 



An updated nearest neighbor search is then conducted: 



30 



35 



40 



45 



50 



55 



and it is ensured that: 



piY\b,A) zKexpi-l^ (y t -b-z t )'-(y t -E-z t )\ 

<l Kexp{-±22 (y c -b-z t )'-(y t -b-z t )\ 
* t 

= p( Y\b, A) . 



Therefore, by iteratively and interleavingly finding the best codeword, z,, and obtaining the best "tentative" 
bias estimate, b, the likelihood function of the bias vector b is increased until a local optimal, or fixed point, 

solution for b is reached. Note that the original, distorted process Y is split into two processes x and B = 

Y - >: . If x is a reasonable estimate of the undistorted signal X, and B is assumed to be stationary, it is 

then reasonable to assume B to be stationary and the maximum likelihood bias estimate b to be a good es- 
timate of the true value of the bias b. Note that an alternative implementation is to gradually reduce the bias 
in the signal at each frame, y t , as the iteration progresses; that is, at the n th iteration: 



5 



BP 0 674 306 A2 



_ „ T 



^ C = l tJ ' 

where z, is defined with; 

* c = y ( - Efcn - i> _ ■ 2).. - 5(0). 
This leads to the same results of maximum likelihood. Also, note that a fixed point solution is reached bfr> -> 
0, as b gets larger. 

For example, a discrete density HMM speech recognition system typically utilizes the generalized Lloyd 
algorithm to compute a locally optimal set of code vectors or centroids to minimize the empirical quantization 
error. In the maximum likelihood formulation of a traditional vector quantizer design, a bias term is included 
which is assumed to be constantly zero. To utilize the SBR method, an identical formulation is used except 
that the centroids are held constant while the bias is treated as an unknown, to be estimated by the maximum 
likelihood method. Therefore, given a set of centroids, u.,, which are computed by the generalized Lloyd algo- 
rithm, the SBR method in a discrete density HMM framework can be carried out for the training and testing 
phases of the speech recognition system. 

Fig. 3 is a flowchart 30 of the implementation of the SBR technique during training of an HMM speech 
recognition system. Atraining speech signal Y(©)). at step 31 , undergoes feature analysis at step 32 to generate 
cepstral coefficients. A first index m is set equal to one at step 33, and a second index n is also set equal to 
one at step 34. An estimate of the spectral bias, b, is computed at step 35 for each utterance of Tf rames, such 
that: 



where the best codeword Zt is the "nearest neighbor" to the distorted signal spectrum y t : 

z t = arsr m J n d(y t/f i t ) . 

One of such distortion criteria consistent with the solution for the bias estimate b is the Euclidean distance: 
Subsequently, the bias estimate b is subtracted from the distorted signal in step 36, so that: 

c = y t -b; 1^mT t 

resulting in a tentative training speech value. * c , which is a maximization of the likelihood function p(xU) 
described above. Next, the index n is checked in step 37 to see if a preset number N has been reached If not 
n is incremented by 1 in step 38, and the bias estimate is recomputed in step 35 using the tentative training 
speech value. This process is repeated until n = N, then step 39 is reached. At step 39 the second index m is 
checked to see if a predetermined number M has been reached. If not, in step 40 the first index n is reset to 
equal one and the second index m is incremented by one. Next, vector quantization is performed in step 41 
and a new codebook having a new set of centroids ( tt ^ is generated in step 42. Once again, a new estimate 

of the bias, b, is computed and then removed from the tentative training speech value , ThuS) thjs entire 

procedure is iterated until M is reached, using the tentative speech value x s rather than y t , to ensure a re- 
duction in the bias (also the quantization error), and a further maximization of the likelihood function. When 
M is reached, the training speech value and optimal centroids obtained are used for training the HMM recog- 
nizer at step 43. In Fig. 3, the dotted line box 45 demarcates those steps which comprise the SBR method of 
the present invention. 

Fig. 4 is a flowchart 50 of the SBR technique used during the testing phase of an HMM recognition system, 
after the HMM recognizer has already been trained. Optimally, the HMM recognizer was trained using the SBR 
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method explained above with respect to Fig. 3. However, other training techniques can be used as long as a 
codebook of centroid vectors was generated. 

Referring to Fig. 4, a contaminated speech signal 51, Y(oo), is input to the system, and feature analysis of 
the signal is performed at step 52 to generate cepstral coefficients. An index n is set equal to one at step 53, 
5 and then an estimate of the bias b is computed as_described above, in step 54, using the codebook 55 that 
was generated during training. The bias estimate b is then subtracted from the speech signal in step 56 to 

generate a tentative speech signal *t . Next, if the index n has not reached a predetermined number N in step 
57, then n is incremented by one in step 58. the bias is recomputed in step 54 and removed from the speech 

10 value to form a new tentative speech signal. This process continues until the predetermined number N Is 
reached at which point the resulting reduced bias signal is fed to an HMM recognizer in step 59 for processing. 
The dottedjine box 46 demarcates the steps of the SBR process. Thus, the equations for computing the bias 
estimate, b, and best codeword, Zt, are repeated with the new, improved set of centroids until the likelihood 
function reaches a fixed point. Typically, one or two iterations are adequate to ensure convergence. Note that 

15 the nearest neighbor search to find the best codeword, A, could instead involve a memory structure such as 
that to be solved by the Viterbi algorithm. 

As discussed above, the distorted speech signal Y(o>) contains two types of biases. These biases can be 
reduced by the above method in an integrated, iterative manner. After minimizing or removing the additive spec- 
tral bias N(a>), the filtered spectral signal may be transformed into cepstrum: 

20 X = x + n. 

where 

* = /DF7Ilog{ ^ (©)}], 
x = /DF7llog{X(©)}] 1 

25 and 

h = /DF7Tlog{H(©)}]. 

By applying the above two-step procedure for bias removal, using cepstrum rather than spectrum, a new 
set of centroids is generated which minimizes or removes the additive bias component This also ensures a 
maximization of the local likelihood probability. The SBR method can then be iterated several times between 
30 the spectral and cepstral domains until extraneous effects are minimized as much as possible. 

It should be noted that there is a strong relationship between signal bias removal and cepstral mean sub- 
traction (CMS). In fact, CMS is equivalent to SBR when a one-codeword vector quantizer is used, where A = 
{uo}. If is a zero vector, thus assuming that the long term cepstral average of speech is zero, then the SBR 
method is reduced to CMS, with the bias vector b representing the frame cepstral average computed over the 
35 whole utterance. 

When incorporating SBR to a platform speech recognizer, an important consideration is the look-ahead 
frame delay necessary for the estimation of the bias. The above discussion assumes that the entire utterance 
is available prior to computing the bias, which is typically not the case for real-world applications. In many prac- 
tical systems, a speech utterance is commonly analyzed on a frame- by-frame basis, or frame synchronously. 

40 A speech frame is equal to some predefined speech interval, for example, 30 millisecond sections of a speech 
utterance. Thus, in real world applications, processing is typically carried out in synchronous fashion wherein 
the first frame is processed by a first microprocessor and then the processed frame passed to one or more 
other microprocessors for recognition analysis as the first microprocessor starts to work on the next frame of 
speech. Thus, acoustic features are passed to the recognizer at every frame, instead of in a batch mode where- 

45 in the entire speech utterance is analyzed all at once. This process of dealing with each frame individually is 
crucial for reat-time implementation and minimal memory requirements. 

Fig. 5A is a block diagram 60, illustrating a speech recognition that can incorporate the present invention. 
A contaminated speech signal Y(oo) is input to a first microprocessor 61 which performs feature analysis using 
software routines stored in a shared memory 62, which comprises both random-access and read-only memory. 

so The first microprocessor 61 also implements the SBR process of the present invention, and speech data Is 

stored in the memory 62. The output speech signal x (©) from the first microprocessor is then input to a sec- 
ond microprocessor 63, which performs speech recognition to generate the text output. 

Fig. 5B is a block diagram 70 illustrating another speech recognition speech recognition apparatus which 
55 may incorporate the invention. A contaminated speech signal Y(co)) is input to a first microprocessor 71, and 
once again a shared memory 72, comprising both random access and read-only memory, is used for data stor- 
age. A plurality of microprocessors (73, 74 to X) process the output speech signal x (to) from the first micro- 
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processor to perform speech recognition resulting in the text output. An apparatus such as that shown in Fig. 
SB is typically used to process a speech utterance on a frame-by-frame basis wherein each microprocessor 
73, 74 to X may operate on different frames as the utterance is processed. Thus, the present invention can 
be utilized as a modular addition to the first microprocessor's routines without the need to modify the process 
used in the speech recognizer and without requiring additional hardware. This modular characteristic is ad- 
vantageous for use with existing speech recognition systems. 

When not all of the frames of the test utterance are simultaneously available for the computation of the 
bias vector b, there are two possibilities to consider. First, a two-pass process may be applied where the bias 
vector as well as other intensive operations are computed in the first pass, leaving the second pass to perform 
recognition. Second, a sequential method for bias removal may be applied such as sequential signal bias re- 
moval (SSBR). If SSBR is to be conducted frame synchronously, then there exist many possible solutions one 
of which is outlined below. 

If bNt. , denotes the average bias vector at the (t -1)* frame, and bt Is the deviation vector at frame tfb. = 
y ( - zj, then: _ * 1 

b t = ab t . , + (1 -a) -b t , 0^a^1, 

and 



Xt = yt - b b 

where a is a weighting coefficient. Note that the estimate of the deviation vector E is computed iterativelv 
2Q such that: Jf 

b t - b t ( n - 1 > + o>- 2 > + .. +Z>/°>, 
where n is the number of iterations for reestimating the bias and b>- D is the sequential bias estimate at iteration 
n - 1 and frame t. In order to ensure that the sequential bias estimate at t = T is equivalent to that computed 
over the whole utterance, the weighting coefficient a is set to (t - 1)/t. Other approaches to performing SSBR 
25 such as using a bootstrapped bias estimate or a leaky integrator are possible. 

It should be realized that the same procedure outlined above with respect to the above equation for the 
deviation vector b t can be used to construct a sequential estimate of the cepstral average. Thus, if q. , denotes 
the average cepstral vector at frame t, and q denotes the cepstral vector at frame t, then: 

c t _ , = a c t . i + (1 - a) Cf, 0^a^1, 

30 and 

t = yt - c t . 

In order to test the invention, a discrete density HMM recognition system was modified to include the SBR 
method (see Fig. 2B). During training of the system, an estimate of the cepstral bias was computed for every 
35 training utterance and subtracted from it. This procedure of estimating and removing the bias was repeated 
twenty times, beyond which no significant reduction in the average norm value, or length, of the bias was ob- 
served. The cepstral codebook was then recomputed and the overall process was iterated four times to ensure 
adequate reduction in the norm of the bias, which was also accompanied by a reduction in the quantization 
error. 

40 The experiments were conducted with an input signal, sampled at 8 kHz, that was initially pre-emphasized 

(1-0.95Z-1) and grouped into frames of 256 samples with a frame shift of 80 samples. Each frame was Ham- 
ming windowed, Fourier transformed into the power spectral domain, and then passed through a set of 30 tri- 
angular band-pass filters. Mel-based cepstral parameters were used, which take advantage of the human au- 
ditory system by sampling the spectrum at mel-scale intervals. In order to compute mel-based cepstral para- 

45 meters, or mel cepstrum, the inverse discrete cosine transform was applied on the smoothed log power spec- 
trum and 12 coefficients were extracted. The first and second order time derivatives of the cepstrum, the delta 
cepstrum and delta-delta cepstrum, were also computed. 

Besides the cepstrai-based features, the log of the energy and its first and second order time derivatives 
were also computed. Thus, each speech frame was represented by a vector of 39 features. Note that the com- 

so putation of all the higher order coefficients was performed over a segment of five windows. 

The input features, namely, 12 cepstrum, 12 delta cepstrum, 12 delta-deita cepstrum, 1 energy, 1 delta 
energy and 1 delta-delta energy were treated separately for Vector Quantization. The generalized Lloyd algo- 
rithm was applied on the entire training data to generate six codebooks, one per feature vector. The codebook 
size was set to 256 for the cepstrum-derived coefficients, and to 32 for the energy-derived coefficients. Such 

55 codebook sizes have been shown to provide a reasonable trade-off between computational complexity and 
recognition performance. The number of iterations for refining the cluster centroids was set to a maximum of 
10, and the generalized Lloyd algorithm employed the L2 norm, with no liftering, for computing the cepstral 
distance. 
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The speech recognizer was based on a discrete density HMM using whole word models, one model per 
digit (1,2..., oh, zero) and per gender, male and female. Models were I eft- to-right with no skip state transitions. 
A total of 24 models including silence and pause were used. The number of states for each model varied be- 
tween one state for silence, and twenty-one states for zero. (The number of states for each digit model was 

5 computed using a simplex search optimization method.) Ten iterations of the maximum likelihood estimation 
were employed during training, followed by three iterations of the maximum mutual information. The latter train- 
ing criterion was also applied for computing the exponents or weights of the six codebooks. The examples as- 
sume unknown length grammar and unenpointed strings both in training and testing. 

Two connected digits databases were used to evaluate the robustness characteristics of the signal bias 

10 removal method. The databases were recorded over telephone lines by having individuals read digit strings 
from a predefined list. 

The first database, DB1, was collected from five dialectically distinct regions within the United States, 
namely, Long Island, Chicago, Boston, Columbus, and Atlanta. Each region consisted of 100 adult talkers (50 
males and 50 females), each speaking 66 connected digit strings from a predefined list (11 digit strings for 
is each of lengths two through seven). Half of their utterances were recorded using two electret microphone hand- 
sets, and the other half using two carbon button microphone handsets. Speech was transmitted over a long- 
distance telephone network that was either all analog, all digital or a mix, depending on the region. A subset 
of this database consisting of 14629 strings was assigned for training, and a different subset of 7073 strings 
was assigned for testing. 

20 The second database, DB2, was collected from two dialectically distinct regions, namely. Long Island and 

Boston, over a digital T-1 interface. Speech was recorded using four different microphone handsets, two elec- 
tret and two carbon button. Digit strings of lengths 10, 14, 15 and 16 digits, corresponding to credit card num- 
bers and long-distance telephone numbers, were collected from 250 adult talkers (125 males and 125 females). 
A subset of this database of 2842 strings was utilized for testing only. 

25 Training was performed on the training portion of the first database D81, and testing was performed on 

the testing portions of both databases DB1 and DB2. Testing on DB1 was considered as under "matched" con- 
ditions, and testing on DB2 was considered as under "mismatched" conditions. 

In order to quantify the degradation in the recognition performance when training and testing in mismatch- 
ed training and testing environmental conditions, the table of Fig. 6 shows the error rate for SBR (column 4) 

30 when using mel cepstrum, for the baseline recognition. The majority of the errors were due to an increase in 
the deletion rate, although a moderate rise in the rates of substitution and insertion was observed. 

Simulation results are described below that illustrate the capabilities of the SBR method, where the com- 
putation of the bias Is performed at the utterance level, and the SSBR method or sequential estimate of the 
bias, when each method is integrated as part of the baseline HMM system. Note that although the formulation 

35 of the bias removal method presented earlier is applicable to the spectral domain for noise bias removal, it 
was strictly employed in the cepstral domain here. 

A plot of the norm of the bias vector | b | at every iteration, averaged over all the training data of DB1, is 
shown in Fig. 7. Note that every time a new code book was generated, the original "unbiased* data was used 
for recomputing the bias. This explains the sudden jump in the norm when m is incremented. When the original 

40 "unbiased 1 ' data or the "biased" data from the proceeding iteration, m - 1 , was used it led to the same results. 
Although using the latter alternative would ensure a faster convergence, it is expensive since additional mem- 
ory storage is then required for the processed data at every iteration. Further, no significant reduction in the 
norm value of the bias beyond m = 2 and n = 10 [n = iteration] was observed (see Fig. 7), and It approaches 
zero as the number of iterations increases. 

45 During recognition, an estimate of the bias was computed for each test utterance and subtracted from the 

speech signal. Similarly this procedure was repeated twenty times. Each utterance was then passed to the 
recognizer. Fig. 7 shows the average norm of the cepstral bias for the test data of DB1 . Clearly, the norm value 
becomes approximately zero beyond ten or so iterations, approaching that of the training data. Fig. 8 shows 
the variations in the quantization error, or the average Euclidean distance, for the training and the testing data 

so at every iteration of the SBR method. The plots indicate a drop in the error by about 30% below its initial starting 
value, before applying SBR, and is only 6% above that of the training data. 

Referring again to Fig. 6, column 4 of the table shows the word error rate when introducing SBR with mel 
cepstrum. These results suggest that the SBR method is able to reduce the word error rate by as much as 
41% for mismatched training and testing conditions (DB2 with mel cepstrum). In addition, the improvement 

55 was chiefly due to a reduction in the rates of deletion and substitution by over a half, during mismatched con- 
ditions, with the insertion and deletion rates becoming relatively equal. 

A similar experiment was also conducted using CMS rather than SBR. Column 6 of Fig. 6 shows the word 
error rates for CMS when using mel cepstrum. The additional improvements that SBR provides over CMS sug- 
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gests that having a finite number of codewords for the computation of the bias that span the acoustic feature 
space, as in the case of SBR, is more rewarding than using a one entry codeword with zero coefficients, as 
in the case of CMS. The relative reductions in the word error rate when introducing SBR, over CMS, for DB1 
and DB2 were 9% and 16%, respectively. A careful examination of the errors that remained following the in- 
troduction of SBR showed them to be a subset of those of the CMS method. 

An advantage of the SBR method over CMS is that SBR can be employed during the testing phase only, 
since the bias estimate, b, is largely dependent on the model, A. Referring to Fig. 6, column 5 (SBR*) of the 
table shows the results of experiments introducing SBR during testing only. As such, the penalty for not training 
with SBRJs rather small. Further, the same approach cannot be applied when using CMS since the cepstral 
average, c. is estimated independently of the model, A. Using CMS during testing only causes an unacceptable 
degradation in the recognition performance (see Fig. 6, column 7, labelled CMS*). 

The effect of SBR on the distribution of the bias is estimated in Fig. 9, which shows a histogram 80, or 
probability distribution, for the second coefficient of the cepstral bias vector when testing on DB1. This dis- 
tribution, as well as those of other coefficients, have been observed to have a Gaussian-like shape. However, 
this was generally not found to be the case for the DB2 database. Fig. 10 shows the distribution 90 for the 
same coefficient when testing on DB2, and depicts a shifted flat distribution with higher variance and dynamic 
range. After several iterations of the SBR method, the bias distributions become sharper (smaller variance) 
and have a smaller dynamic range. This is illustrated in Fig. 11, where the histogram 95 maps to that of Fig. 
10 following ten iterations of the SBR method. 

All of the previous experiments assumed the existence of the entire utterance for the estimation of the 
bias, therefore the look-ahead frame delay was equivalent to the size of the input utterance. 

As explained above, to realize the SBR method in a frame synchronous mode, a sequential estimate of 
the bias (SSBR) can be applied for recognition. The SSBR method includes an on-line estimate of the bias 
that is regularly updated on every frame using the criterion set forth in the equation for the b t bias vector. In 
the examples which follow, ten iterations were performed for each frame for the computation of the bias vector, 
b t) prior to updating the input signal, *c 

Fig. 12 presents the word error rates when using mel cepstrum for the baseline system (column 2) with 
eitherSBR (column 3), orSSBR (column 4), or CMS (column 5). Note that the SSBR method was strictly applied 
during testing, while training included the application of the SBR method (i.e., results for both the SBR method 
and the SSBR method utilized the same training models, but differed in the manner of computing the bias dur- 
ing testing). Clearly, the penalty for having a sequential estimate of the bias, as opposed to estimating the 
bias in a batch mode, is rather small. The relative changes in the word error rate when introducing SSBR, over 
SBR, for testing on DB1 and DB2 were +1 5% and +6%, respectively. 

One attractive feature of the SBR method is that its formulation does not contain parameters that have 
to be tuned differently for different databases. The SBR method indicates that by iteratively estimating and 
subtracting the bias for each individual utterance in the manner presented above, a maximization of the like- 
lihood probability function results. 

Another interesting feature of the SBR method when integrated as part of a VQ (or a semi-continuous)- 
based speech recognition system is that one can apply the same cepstral codebook used by the HMMs for 
the computation of the bias, thus requiring no additional memory storage. 

Fig. 1 3 shows the word error rate for the SBR method as a function of the codebook size, where the co- 
debook size for the computation of the bias was varied from 1 to 1 024 entries when testing on DB1 and DB2. 
These results indicate that a codebook size of 3 to 4 bits (or 8-16 which equals 2*~ 1 to 2^) may be sufficient 
for the computation of the bias without significant loss in performance. 

Fig. 14 shows the word error rate as a function of the codebook size when employing the SSBR method. 
Introducing a small-sized codebook in this case has a disastrous impact on the recognition performance since 
the noisy estimate of the bias during the initial part of the utterance causes removal of useful speech informa- 
tion, rather than extraneous signals. Therefore, a much larger codebook (e.g., 128 to 256) is necessary for 
the SSBR method to have a positive impact on the recognition performance. Thus, when applying sequential 
signal bias removal to the speech signal, a moderate to a large sized codebook should be used if a minimal 
look-ahead frame delay is desired. If look-ahead frame delay is of no concern then one may use a small sized 
codebook (this is the case for SBR). 

A related problem is the issue of what size string is necessary for "adequate" computation of the bias. Fig. 
15 is a table which shows the string error rate before and after the application of the SBR method when using 
mel cepstrum. Although the testing is performed on DB1 under matched conditions, it suggests that a minimum 
of a three digit string (about 1 second) is necessary in order for the bias removal method to have any positive 
impact on the recognition performance. Experiments conducted on single digit strings show that the current 
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formulation of the SBR causes a slight degradation of up to 5% in the digit error rate. 

In summary, the signal bias removal (SBR) method utilizes an iterative procedure for estimating the bias 
in the spectral and cepstral domains for the minimization of deleterious signal components in telephone speech 
recognition. The procedure is based on maximizing the likelihood of a speech mode! in which the bias is con- 

5 sidered as the unknown parameter. The SBR method, as applied in the cepstral domain only, can be integrated 
as part of a discrete density HMM system. Further, to enable real-time implementation, a sequential signal 
bias removal method (SSBR) was shown to be effective when processing speech signals on a frame-by-frame 
basis. Results from experiments using two speaker-independent databases, wherein the data from the speak- 
ers consisted of spoken strings of digits, indicate that the SBR method, when applied to a fairly long string of 

w digits, is capable of minimizing extraneous channel distortion, and consequently improving the performance 
of telephone speech recognition. 

Further, the experimental results indicate that when introducing SBR during testing only, as opposed to 
during both training and testing, the word error rate only rises up to 14%. For CMS, this would result in a jump 
in the error rate by a factor exceeding three times. This advantage of being able to apply SBR without retraining 

15 the recognition models is desirable in all existing applications of speech recognition. 

It is to be understood that the above-described embodiments are merely illustrative, and that many vari- 
ations can be devised by those skilled in the art without departing from the scope of the invention. 



20 Claims 

1. A method for minimizing the effect of an unknown signal bias in an input speech signal for use by a speech 
recognition system, comprising: 

(1) training the speech recognition system by using the following steps: 
25 (a) generating a set of centroids based on a training speech signal; 

(b) computing an estimate of the bias for the training speech signal based on maximizing a likelihood 
function; 

(c) subtracting the estimate of the bias from the training speech signal to obtain a tentative training 
speech value; 

30 (d) repeating steps (b) and (c) a preset number of times, wherein each subsequent computed esti- 

mate of the bias is based on the previous tentative training speech value to arrive at a reduced bias 
training speech signal value; 

(e) recomputing the centroids based on the reduced bias training speech signal to generate a new 
set of centroids; 

35 (f) repeating steps (b) to (e) a predetermined number of times to compute a processed reduced bias 

speech signal and to form an optimal set of centroids; 

(g) utilizing the optimal set of centroids and the processed reduced bias speech signal as training 
input for a speech recognizer; 

(2) testing an input speech signal to minimize the unknown bias by using the following steps: 

40 (h) utilizing the optimal set of centroids to compute an estimate of the bias for each utterance of 

the speech signal based on maximizing a likelihood function; 

(i) subtracting the estimate of the bias from the speech signal to obtain a tentative speech value; 
G) repeating steps (h) and (i) a preset number of times, wherein each subsequent computed esti- 
mate of the bias is based on the previous tentative speech value, resulting in a reduced bias speech 
45 signal value; and 

(3) utilizing the reduced bias speech signal as input to a speech recognizer. 

2. The method of claim 1 , wherein the speech recognition system utilizes a Hidden Markov Model speech 
recognizer. 

50 

3. A method for minimizing the effect of an unknown signal bias on an input speech signal during the testing 
phase of a speech recognition system, comprising: 

(a) computing an estimate of the bias for each utterance of the speech signal based on maximizing a 
likelihood function by initially utilizing a set of centroids generated by a training model; 
55 (b) subtracting the estimate of the bias from the input speech signal to obtain a tentative speech value; 

(c) repeating steps (a) and (b) a predetermined number of times, wherein each subsequent computed 
estimate of the bias is based on the previous tentative speech value, resulting in a reduced bias speech 
signal value; and 
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(d) utilizing the reduced bias speech signal value as input to a speech recognizer. 
4. The method of claim 3, wherein a vector quantization method is utilized to generate the centroids of step 



5. 



A method for sequentially minimizing the effect of an unknown signal bias of an input speech signal for a 
speech recognition system, comprising: y 
(1) training the speech recognition system by using the following steps: 

(a) generating a set of centroids based on a training speech signal; 

(b) analyzing the speech signal on a frame by frame basis or in a batch mode; 

(c) computing an estimate of the bias for the training speech signal based on maximizing a likelihood 
function; 

(d) subtracting the estimate of the bias from the training speech signal to obtain a tentative training 
speech value; a 

(e) repeating steps (c) and (d) a predetermined number of times, wherein each subsequent com- 
puted estimate of the bias is based on the previous tentative training speech value to arrive at a 
reduced bias training speech signal value; 

(f) recomputing the centroids based on the reduced biased training speech signal value to generate 
a new set of centroids; y 

(g) repeating steps (c) to <f) a predetermined amount of times to compute a processed reduced bias 
speech signal and to generate an optimal set of centroids; 

(h) utilizing the optimal set of centroids and the processed reduced bias speech signal as training 
input for a speech recognizer; 

(2) testing an input speech signal to minimize the unknown bias by using the following steps- 

(i) setting a weighting coefficient for updating a bias value; 
(j) analyzing an utterance on a frame- by-frame basis; 

(k) computing a sequential bias estimate for each frame of the speech signal based on maximizing 
a likelihood function; a 

(I) subtracting the sequential bias estimate from the input speech signal at every frame to obtain 
a tentative speech value; 

(m) repeating steps (k) and (1) a predetermined number of times, wherein each subsequent com- 
puted estimate of the bias is based on the previous tentative speech value, resulting in a reduced 
bias speech signal value; and 

(3) utilizing the reduced bias speech signal as input to a speech recognizer. 

. The method of claim 5 f wherein the speech recognition system utilizes a Hidden Markov Model speech 
re cog n izer. 

. A rnethod for sequentially minimizing the effect of an unknown signal bias on an input speech signal during 
the testing phase of a speech recognition system, comprising: 

(a) setting a weighting coefficient for updating a bias value; 

(b) analyzing an utterance on a frame-by-frame basis; 

r? ^TT 9 8 seq " entia ' WaS eStimate f0r 83011 frame of the s P eech si 9 naI base d on maximizing a 
likelihood function by utilizing a set of centroids generated by a training model- 
ed) subtracting the sequential bias estimate from the input speech signal at every frame to obtain a 
tentative speech value; 

(e) repeating steps (c) and (d) a predetermined number of times, wherein each subsequent computed 
estimate of the bias is based on the previous tentative speech vaiue, resulting in a reduced bias speech 
signal value; and 

(0 utilizing the reduced bias speech signal as input to a speech recognizer. 
The method of claim 7, wherein a vector quantization method is utilized to compute the centroids of step 

A speech recognition apparatus for minimizing the effect of an unknown signal bias on an input SDeech 
signal, comprising: K H 

a first microprocessor means for performing feature analysis on the input speech signal and for 
performing a signal bias removal technique to generate a reduced bias signal; 
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a shared memory means connected to the first microprocessor for storing speech data; and 
a second microprocessor means connected to the memory and connected to the output of the first 
microprocessor means, for performing speech recognition based on the generated reduced bias signal. 

The apparatus of claim 9, further comprising: 

a plurality of microprocessor means connected to the output of the first microprocessor and con- 
nected to the memory, for collectively performing speech recognition based on the generated reduced 
bias signal. 
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